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Abstract 

We speculate that the universe is filled with stars composed of electromag- 
netic and dilaton fields which are the sources of the powerful gamma-ray bursts 
impinging upon us from all directions of the universe. We calculate soliton- 
like solutions of these fields and show that their energy can be converted into 
a relativistic plasma in an explosive way. As in classical detonation theory 
the conversion proceeds by a relativistic self-similar solution for a spherical 
detonation wave which extracts the energy from the scalar field via a plasma 
in the wave front in the atmosphere of the star. 



1 Introduction 

Many modern theories of the universe assume the existence of various types of scalar 
fields. Such fields could explain the recently discovered acceleration of the expansion 
of the universe (see e.g. Ref. |TJ), or the formation of clustered systems leading to 
*email: astra@freenet.kg Physics Institute of NAN KR, 265 a, Chui str., Bishkek, 720071, 
Kyrgyz Republic 

^email: kleinert@physik.fu-berlin.de Inst. f. Thcor. Physik, Freie Univ. Berlin 

■femail: hjschmi@rz.uni-potsdam.de Inst. f. Math. Univ. Potsdam, Am Neuen Palais 10, 

D- 14469 Potsdam, Germany 



1 



gravitational walls for galaxies and galaxy clusters. If such fields really exist, the 
universe could contain many compact star-like configurations of large total mass, 
called scalar stars. 

Among the many possible scalar fields, the dilaton field has a special theoreti- 
cal appeal. It couples in a unique minimal way to electromagnetism to make the 
Maxwell action dimensionless. This coupling leads to star-like objects which are 
composed of scalar and electromagnetic fields. That such objects can exist was 
pointed out by many authors [@, ||, |J. We argue that such electro-dilaton stars may 
be responsible for the strong gamma-ray bursts observed in the universe. 

At present there exists no simple conventional explanation for the origin of these 
events reaching us isotropically from all directions of the universe. For an compre- 
hensive discussion of the subject, in particular for the failures of most theories, see 
the article by Ruffini ||. At the same time, there exist models for the early uni- 
verse where an initial state of large volume of relativistic plasma quickly expands 
as a spherical wave, like a critical bubble in an overheated liquid, causing a decay 
of a false vacuum, and creating the universe from such a bubble ||. The electro- 
dilaton stars can supply such an explanation. We assume that the dark matter in 
the universe contains a multitude of such objects, whose total mass exceeds by far 
the total mass of luminous matter and is responsible for the large-scale structure of 
the universe. The luminous matter concentrates in the gravitational potential wells 
of the scalar fields. 

Let us imagine that collision of relativistic particles produce a fireball of critical 
size. Such bubbles have been investigated as possible triggers of phase transitions 
in the early universe ||. In the context of gamma-ray bursts, similar assumptions 
have been made in Ref. 0. The critical bubble forms the seed for transferring 
effectively scalar fields into pairs of elementary particles and their antiparticles. The 
transfer may be initiated by fast oscillation of a field on the outer boundary of the 
fireball. The process causes a relativistic detonation. In a conventional detonation, 
chemical energy is converted into kinetic energy. In a relativistic detonation of 
an electro-dilaton star, it is the energy of the electric and scalar fields which is 
rapidly converted in particle-antiparticle pairs. The resulting fireball expands with 
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relativistic velocity and will therefore not depend on the weak gravitational fields of 
a Newtonian configuration. We may thus study the process within special relativity. 



2 Relativistic Detonation 



We begin by deducing the self-similar solutions for a spherical relativistic detonation 
which goes over to the well-known Zeldovich solution in the limit of small veloci- 
ties ||. The set of equations of relativistic hydrodynamics is conveniently described 
in a spherical coordinate system r, G, <fi. If v denotes the radial velocity of the plasma 
in three-dimensions || [| and e is energy density, p - pressure, the equation of motion 
reads 
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where 7 2 = 1 — v 2 , W = e + p, and c = 1. As in the nonrelativistic case, the motion 
of the plasma behind the detonation front is considered as isentropic, such that ([I]) 
and (0) are the only relevant equations. 

The pairs of relativistic particles and antiparticles created behind the wave front 
generate a high-temperature plasma with the equation of state: 
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where c s is the sound velocity. As in an ordinary spherical detonation problem, we 
search for a solution depending on the self-similar variable 
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in which the differential equations (ED) and (0) reduce to ordinary differential equa- 
tions, which can be combined to the equation for v: 
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In the nonrelativistic limit where v(£), £ 1, this equation reduces to a textbook 
equation M. The qualitative analysis of the relativistic equation (M is similar to 



the one in the textbook: The solutions for v and e have a diverging derivative at 
the wave front. Such singularity of a derivative is defined by the following: the 
expression in parentheses is velocity of current of plasma in relation to the wave 
front. According to the theory of detonation, this velocity is equal to sound velocity 
c s . Therefore, at approaching to the wave front the expression in brackets tends to 
zero. The coordinate of the wave front £ = (where Vd is a velocity of detonation 
wave in central frame) and the velocity of gas behind the wave in laboratory frame 
v g is uniquely determined from conservation laws at the wave front. 
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Figure 1: Dependence of the energy density e and velocity v of plasma on the 
self-similar variable £ behind the front of detonation wave. 



The discussion of possible mechanisms of "recycling" of the field behind the wave 
front into relativistic plasma is considered below. Let us specify here the following 
estimates of the value of and energy density behind the detonation wave. Consider 
a scalar field in the simplest form, with an energy-momentum tensor 

1 



-<p tli (p>» - V(<p) 



V(<p) = mV/2. 



(6) 



In the "scalaron" regime where it undergoes fast oscillations of frequency m, i.e. 
where ip{r, t) = a(r) smmt, the spatial gradients of the field can be ignored and the 
energy density in the laboratory frame) is 



£ f « m 2 a /2. 



(7) 
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The expression for Vd and the energy density of the plasma behind the wave front 
are determined from the conservation laws Tq (field) = Tq 1 (plasma) and T± (field) = 
T± (plasma) for an observer moving with the wave front. Recall that the plasma is 
emitted from the wave front with the velocity c s . Hence 

2c, 2 
v g = c s - v d = e p = Y^ £ f- ( 8 ) 

Since a relativistic plasma has c s = l/v3, we find Vd = v3/2 and £p — 3 £/. The 
self-similar solutions for this case are presented in Fig. |I[ 

We now turn to the mechanism of transition of the field energy into a relativistic 
plasma. 



3 Electro-Dilaton Wave 

So far, the description of the conversion of the energy of a scalar field into relativistic 
plasma at the front of "detonation" wave is purely phenomenological. The physical 
properties of the front of the detonation wave are completely determined by the 
conservation laws for a general energy-momentum tensor Tf . 

One specific mechanism of such a conversion was considered in |12|] based on an 



analogy with a laser. Here we shall consider an alternative possibility where the 
relativistic plasma and radiation are obtained from the energy of a dilaton field. 
The star-like configurations for such field with strong and weak (Newtonian) 



gravitational field have been considered before [TJ], [15|, [HJ. As in Section [| we 
shall discuss only the case of a weak gravitational field and treat the problem within 
special relativity. 



The Lagrangian density of a system of dilaton and electromagnetic fields is [13 



L = 2<S>, i <S» i -((F lm F lm )e- 2a *. (9) 

The parameter ( can have the values ±l[]as will be explained below. The normaliza- 
tion of the fields is the same as in [0, [15[ . The unique interaction of the scalar field 
$ with an electromagnetic field required by scale invariance allows for a nontrivial 
1 our choice is £ = — 1 
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combined electro-dilaton configuration. In the front of the detonation wave, the 
electromagnetic field reduces the dilaton field strength by dissipation. Depending 
on the intensity of the electric field, the dissipation may be due to the creation of 
pairs of particles and antiparticles and to a heating of the plasma. The combined 
process — generation of an electromagnetic field and its subsequent dissipation — 
supplies the energy to the front of the detonation wave. Let us study the process 
in a simple plane-wave configuration. The equations for $ and F^. following from 
Lagrangian (||) are 



= 0, (10) 
#1 = -^-e- 2a *(F lm F lm ). (11) 

The system is supplemented by the missing Maxwell equations (electromagnetic 
versions of the Bianchi identities): 

e Mm F u , m = 0. (12) 
The total energy-momentum tensor associated with the Lagrange density @ is 

T k = 2 $ _ 2(e- 2a ' s 'F a F lk - - #[2($ ,$•') - (e~ 2a * (F lm F lm )]. (13) 



4 Longitudinal Electro-Dilaton Wave 

We now show that there exist plane electro-dilaton waves travelling along the x-axis 
in which the electric field has a component F 10 = E x with all other components 
vanishing. Whereas Eq. (|i~2"|) is fulfilled trivially, Eq. (0) yields 

(e- 2a *E x ), T = 0; (e- 2a *E x ) >x = 0. (14) 

Thus we find a constant of motion 

e- 2a ^E x = E = const. (15) 

The quadratic field combination 

I = F lm F lm = -2E 2 X (16) 
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has the same negative value in all systems of special relativity. For this reason we 
follow Ref. [14] in choosing the parameter ( = — 1 in (H). 
The other field equation ( |TTD becomes 

<9$ <9$ 



aE 2 e 2a ®. 
Or 2 dx 2 



(17) 



or, introducing new variables 2aE x ~~ > 2aE r — > r and function ip = 2a$, we 
have 

^ ^ e*/2. (1> 



<9r 2 <9x 2 



This is one of form of Liouville equation [0. Its possesses a steady-state wave 
solution which, after the redefinition of the variables 



£ = x — ur, 77 = £/ vl — u 2 



takes the form 



c/ 2 </> 
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The first integral of this equation leads to the differential equation 

2 
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whose solution 



(19) 



(20) 



; 2 i) 



-21ncosh(r//2), £ x = E / cosh 2 (r//2) 



(22) 



describes a soliton. Though the potential of dilaton field ip diverges at \r)\ — > 00, 
where $ ~ — 77, the derivatives $ t j remain finite, thus ensuring a finite energy density 
of the field for all rj. Note that in contrast to a charged plane in electrostatics, the 
electric field of this solution has a zero flux E x at infinity. Asymptotically, no electric 
field is detectable. 

There exists a general class of solutions of the Liouville equation fll8| ) containing 
two arbitrary functions fi(x — r), ^(x+r). Setting ip* = 2aQ + \iaaE 2 , the solution 
of eq. (0) is 



ip* = In 



cosh 2 [/i(x - r) + f 2 {x + r)] 



(23) 
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The primes denote derivatives. Especially simple solutions are obtained for the 
linear functions 

1 + u 

fi{x-r) = j(x-r), 7= _ , (24) 

1 — 71 

/ 2 (x + r) = /?(x + r), /?= - (25) 

4y 1 — m z 

where u is velocity of soliton in Eg. (|i"9|) . 

The following solutions are of special interest: 

a) Localized solution for E x with a fixed asymptotic energy density e(|x| — > 
oo) = e 0) where goes over into the previous soliton solution ( p2]) . In these 
solutions, the functions /x and /2 in (23) and their derivatives are regular. 



Singular solutions of Liouville equation are also known [18|. They may lead 



to a local catastrophic growth of electric field, and require special attention. 

b ) The main purpose of this section is to show that at different values of the 
energy of electro-dilaton wave before and behind the electric layer the electro- 
magnetic energy may increase in time. A similar distributions of the energy 
may be created by gravitational configuration. 

Let us specify the arbitrary functions in Liouville solution as 

fi = 7(x-r)+F, F = -lntanh/x(a; - r), (26) 
h = P(x + t), (27) 



where 7, {3 are determined by the initial velocity of the wave from (|24j ) and 
([25]), and \i is an arbitrary constant. Now the profile of the electric field looks 
like 

Ex = 16/3[7 + /jtanh//(r-a;)] 

E cosh 2 [(7 + (3)x — (7 — (5)t + In coshyu(r — x)] 

The profile of such a wave is shown on Fig. 0. Let us assume that the numerator 
of Eq. (123) slowly changes in comparison with the denominator. In that case, 
at given r the maximum of electric soliton concentrates about a zero point of 
argument of cosh. As we see from Eq. (|24j), this region moves with velocity 
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Figure 2: Growth of electrical soliton in wave with different asymptotic values of 
energy both before and behind wave. Solutions (|28|) with u = l/\/3 and [i = 7 
are shown for increasing times r. The energy density in units of E 2 is equal to 7.4 
behind and 0.2 in front of the wave front. For large r the solution gives stationary 
soliton with smaller value of energy. 

u < 1, so that the argument of tanh is always positive in this region. The 
numerator and thus the electric field increase with time. This increase comes 
to an end as tanh[/i(r — x)} reaches its asymptotical value 1. This is illustrated 
in Fig. |2|. If instead of tanh[/i(r — x)] in (|26|) a function is chosen which grows 
without bounds, then also the electric field will keep growing. 

The growth of the energy of the electric field will be consumed by dissipation. 
Its influence on the electro-dilaton wave will be considered in Section |[ 



5 Wave of Transverse- Magnetic Type 

Let us use the above results to study a wave with a transverse magnetic wave. In 
this case the Maxwell equations ([12]) are not fulfilled identically but yield a relation 
between transverse components F^. Let us select for consideration the following 
nonzero components F 20 = E y ,F 21 = H z . Then the complete set of the equations 
is 

Ue-^E y ) + JV 2 ^) = 0, (29) 
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dx dr 
A travelling wave has the field components 



dE v dH z , . 

" + ur = °- (3°) 



E y = uH z (C); H z = H,e 2a "- I = 2H 2 (1 - u 2 )e 4a " . (31) 



Since u < 1 and 7 > 0, the field is transverse magnetic. It dictates in (10) the choice 



of the sign £ = 1. Further, using similar variables as in ( fL9|) 

f = x - ut, ip = 2a$, 77 = 2a£\H \, (32) 
we obtain the solutions for dilaton and electromagnetic fields similar to (|22| ) 

iff = -21ncosh(7?/2); E y = uH z = uH / cosh 2 [r]/ 2). (33) 
The principal difference with respect to (p2|) is the absence of a relativistic factor 



vl — v? in the argument 77. 

The transverse magnetic wave is focused in a band of width Ax ~ 1/o\Hq\. 
Outside of this, the energy density lies mainly in the gradient part of the dilaton 
field. From Eq. ( |3"3"D it follows that energy flux density is independent of £: 

T° = H 2 (l + u 2 ); n = 2H 2 u. (34) 

This can be interpreted as follows: the conversion of the energy of the dilaton field 
from the gradient to the potential part implies a conversion of the gradient energy 
of the dilaton field to the energy of electromagnetic field. Certainly, this conversion 
is completely reversible. 



6 Dissipation 

The above electro-dilaton soliton appears at the front of the detonation wave. If we 
use the plasma behind the front as a frame of reference, then this soliton will be 
move with a velocity u — c s — 1/ a/3 in the positive x-direction. Such a movement 
generates heat in the plasma and creates of pairs of particles and antiparticles in 
high concentration. This reduces the E and H fields. In the simplest description, 
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the decay wave can be taken into account in a travelling-wave approximation by 
replacing Eqs. ( p9|) by 

(1 - u 2 )^-(e- 2a *H z ) = fe~ 2a *H z ; {e~ 2a *H z ) = # ex P [/£/(l - u 2 )\, (35) 

where / is the friction constant with the dimension of a reverse length. While moving 
though the plasma with £ < 0, the electromagnetic wave decays. The equation for 
the dilaton field with £ = 1 becomes 

Substitution here (|35"D and using the redefinitions (0), we obtain 

H z = H h( V )e^, with (3 = f/2aH , (37) 
^ = -\h\^ = -\e^. (38) 

After a change of the variable if) + 2/?r^ — > if), this equation reduces to the Liouville 



form fl2ip. Then subject to dissipation we have the following solution: 

</> d = -2\j3rj + lncosh^/2)]; H d {rf} = / cosh 2 ( V /2). (39) 

Here, the subscript d denotes a solution with dissipation. Thus electro-dilaton soli- 
tons with dissipation become asymmetric with a steeper front part. The creation of 
particles and heat of plasma happens at decreasing of the energy density of dilaton 
field. For \rj\ —>■ oo the energy of the system is concentrated only in the gradient 
part of dilaton field (see I®,©) 



T« = {l + u 2 )Hl{if, din y. 

In the limit \q\ — > oo, the equations (|39| ) becomes 

T °(oo) = (1 + u 2 )H 2 Q (l + 2f3) 2 ; T °(-oo) = (1 + u 2 )H 2 (l - 2(5) 2 . 

It is clear from here that the energy density of dilaton field before the wave is more 
than the energy density behind the wave (J3 > 0). The limiting value for j3 in 
this example is f3 — > 0.5. This implies that the whole energy of the dilaton field 
transforms to the heat energy of plasma. 
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7 Induced Current 



Consider now the dissipation of energy by the induced fields at the wave front. For 
this purpose we substitute a current on the right hand side of Eq. fl29|) 

3y = ° E y = ™H z (£), 

where a is specific conductivity of the medium. Then we find the following system 
of dimensionless equations: 

with the parameter (3 ~ au. At /3 = 0, this reduces to the previous electro-dilaton 
soliton (|33f) . 

To analyze Eqs. (|40|) it is convenient to introduce the new variable 



z = /ie"^, (41) 

which is unity in the absence of dissipation. The first integral of the differential 
equations (EO) leads to the solution: 



The integration constant is selected so that in the limit r] — > — oo the solution tends 
to the dissipation- free soliton with (dip / drj) _ — > 1 for — > 1. In the opposite limit 
7] — > oo, the variable z 2 tends to zero and 

Thus for z 2 — > 0, the solution represents a kink moving in the positive x-direction. 
By virtue of (3 > 0, the asymptotic value of the dilaton energy density before the 
front is less than behind it. Physically this means that the pressure of the dilaton 
field creates a plasma at the wave front pushing it ahead. It follows from Eq. (p2|) 
that the limiting value of f3 is now (3 C = 1/4 implying that the electro-dilaton energy 
goes completely over into an energy flux of moving plasma. The solutions are plotted 
in Fig. |. 
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Figure 3: Dependence of the energy density e and velocity v of plasma on the 
self-similar variable £ behind the front of detonation wave for various parameters (3. 

8 Conclusion 

^From our discussion it appears perfectly plausible that dark matter consists of 
electro-dilaton stars and is not dark at all, but has been showing its presence quite 
dramatically all along via the powerful gamma-ray bursts observed since. 



9 Acknowledgement 

VF and VG than the DAAD in Bonn/Germany for financial support. We also thank 
Dr. A. Pelster for discussions. 



References 



[1] V. Sahni and A. Starobinsky, ^r-qc/ 990439^ , Int. J. Mod. Phys. D 9 (2000) 



373; a review for the older results is: S. Gottlober, V. Miiller, H.-J. Schmidt, 
A. A. Starobinsky, Models of chaotic Inflation, Int. J. Mod. Phys. D 1 (1992) 
257-279. 



13 



See also the related discussion of quintessence in R.R. Caldwell, R. Dave, P. J. 
Steinhard, Cosmological Imprint of an Energy Component with General Equa- 
tion of State, Phys. Rev. Lett. 80, 1582 (1998) ( |astro-ph/ 9708069) ); 
For experimental data see 

P.M. Garnavich et al., Supernova Limits on the Cosmic Equation of State, As- 
trophys. J. 509, 74 (1998) ( |astro-pli/9806396| ); 

B. Ratra, P.J.E. Peebles, Cosmological Consequences of a Rolling Homogeneous 
Scalar Field, Phys. Rev. D 37, 3406 (1988); 

C. Kolda, W. Lahneman, Exponential Quintessence and the End of Ac- 
celeration, |hep-ph/0105300|. M. Gasperini, Dilatonic Interpretation of the 



Quintessence? , |gr-qc/0105082| ; 

J. M. Cline, "Quintessence, Cosmological Horizons, and Self- Tuning" , [hep- 
ph/0105251 



Ph. Brax, A. C. Davis, Cosmological Evolution on Self-Tuned Branes and the 
Cosmological Constant, JHEP 0105, 007 (2001) (|hep-th/0104"023|) ; 
A. Riazuclo and J. P. Uzan, |astro-ph / 1 07386} 

Bin Chen and Feng-Li Lin, Holographic Cosmic Quintessence on Dilatonic 
Brane World, [hep-ph / 1 60 54| ; M. Gasperini, F. Piazza, and G. Veneziano, 
Quintessence as a run- away dilaton, |gr-qc/0108016 . 



[2] R. Ruffini and S. Bonazolla, Phys. Rev. 187, 1767, (1969). 
[3] T. Piran, Phys. Rev. D36, 3575, (1987). 

[4] E. Mielke and F. Schunck, kr-qc/9801063t Gen. Rel. Grav. 33, 805, (2001). 



[5] R.J. Ruffini, Analogies, New Paradigms, and Observational Data as Growing 
Factors of Relativistic Astrophysics, in Fluctuating Paths and Fields - Dedicated 
to Hagen Kleinert on the Occasion of his 60th Birthday, Eds. W. Janke, A. Pel- 
ster, H.-J. Schmidt, and M. Bachmann (World Scientific, Singapore, 2001), p. 
771. 

See also 

R. Ruffini, C.L. Bianco, P. Chardonnet, F. Fraschetti, S.-S. Xue, On the inter- 
pretation of the burst structure of GRBs, Astrophys. J. 555 (2001), L107, L117; 

14 



Nuovo Cimento B, 116, 99; 

According to these works, gamma-ray bursts are correlated with supernovas. 
The discussion in the present paper attempts to give an alternative interpreta- 
tion. 

G. Mao, S. Chiba, W. Greiner, K. Oyamatsu, Vacuum discharge as a possi- 
ble source of gamma-ray bursts, Int. J. Mod. Phys. E 9, 185 (2000) ( |astro- 
ph/9912303|) . 



[6] S. Coleman and F. De Luccia, Phys. Rev. D21, 3305, (1980). 
[7] K. Postnov, Physics Uspekhi 42, 469, (1999). 

[8] L. Landau and E. Lifshitz, Hydrodynamics (Nauka, Moscow, 1988), P. 680. (In 
Russian) 

[9] F. Baum, S. Kaplan and K. Staniukovich, Introduction in Cosmic Gas Dynam- 
ics (Moscow, 1958). 

[10] A. Linde, Elementary Particle Physics and Inflationary Cosmology (Nauka, 
Moscow, 1990). 

[11] A. Migdal and V. Krainov, Approximate methods of quantum mechanics 
(Nauka, Moscow, 1966). 

[12] I. Tkachev, Phys. Lett. B191, 41, (1987). 

[13] Zhi-jian Tao and Xun Xue, Phys. Rev. D45, 1878, (1992). 

[14] R. Gregory and J. Harvey, Phys. Rev. D47, 2411, (1993). 

[15] T. Matos, D. Nunez and H. Quevedo Phys. Rev. D51, R310, (1995). 



[16] S. Yazadjiev, P. Fiziev, T. Boyadjiev and M. Todorov, [hep-th/0105165 . 
[17] J. Liouville, J. math, pures et appl. 18, 71, (1853). 

[18] G. Jorjadze, A. Pogrebkov and M. Polivanov, ICTP prep. IC/78/126. (Trieste, 
1978). 



15 



